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Introduction 

The control of chaos represents a very real and impor- 
tant problem in a wide variety of applications, ranging 
from neuron assemblies to lasers and hydrodynamic sys- 
tems (jy. The procedure used consists of stabilizing an 
unstable periodic orbit by the application of precisely de- 
signed small perturbations to a parameter and/or a tra- 
jectory of the chaotic system. Different methods of chaos 
control have been suggested and applied in many differ- 
ent physical contexts, as well as numerically to model 
systems 0. For practical applications of these control 
methods, it is important to understand how noise influ- 
ences the stabilization process, because fluctuations are 
inherent and inevitably present in dissipative systems. 
The problem has not been well studied. Typically, a 
method is developed for stabilization of the orbit with- 
out initially taking any account of fluctuations. Only 
then do the authors check the robustness of their method 
by introducing weak noise into the system P|. Thus, in 
the celebrated pioneering work of Ott, Grebogi and York, 
"Controlling chaos" 0, the authors just noted that noise 
can induce failures of stabilization. 

In several works 0, 0] methods are developed for the 
stabilization of unstable orbits in the presence of noise. 
They are based on a strong feedback approach to sup- 
press any deviation from the stabilized states. There are 
also methods 5j that use noise to move the system to a 
desired unstable state, and then stabilize it there. 

In this work we consider noise-induced failures in the 
stabilization of an unstable orbit and the problem of 
controlling these failures. The method of Ott, Grebogi 
and Yorke (OGY) and a modification of the adaptive 
method (ADP) [l| are used to stabilize an unstable point 
of the logistic map. We consider the small noise limit 
where stabilization failures are very rare and can there- 
fore they be considered as large fluctuations (deviations) 
from a stable state. We study the properties of large 
deviations by determining the optimal paths and the op- 
timal fluctuational forces corresponding to the failures. 
We employ two methods to determine the optimal paths 
and forces. The first of these builds and analyzes the pre- 
history probability distribution to determine the optimal 
path and optimal force |3| . The second method considers 



an extended map (relative to the initial one) which de- 
fines fluctuational paths and forces in the zero-noise limit 

Q . Furthermore we use the optimal paths and forces 
to develop methods of controlling the large deviations, 
i.e. the noise-induced failures of stabilizationp^ . 

In section I we describe the procedures for local and 
global stabilization of an unstable orbit of the logistic 
map. The general approach to the control of a large de- 
viation is presented in section II. Noise-induced failures 
of local and global stabilization are considered in sec- 
tions III and IV respectively. The results obtained are 
discussed in the conclusion. 



I. CHAOS STABILIZATION 

For simplicity we will stabilize an unstable fixed point 
X* of the logistic map: 



(1) 



where a;„ is a coordinate, n is discrete time and r is the 
control parameter that determines different regimes of 
the map's behavior The coordinate of the fixed point 
X* is defined by the condition: Xn+i — Xn, and conse- 
quently its location depends on the parameter r: 
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a;* = 1 - -. 

r 



(2) 



We set the parameter r — 3.8, a value for which an ape- 
riodic (chaotic) regime is observed 1^, and the point x* 
is embedded in the chaotic attractor. 

From the range of existing stabilization methods, we 
chose to work with just two: the OGY and ADP methods 
mentioned above. 

To stabilize a fixed point by the OGY method, pertur- 
bations Ar are applied to the parameter r, leading to the 
map being modified Q in the following manner: 

Xn+l = {r + Arn)Xn{l~ Xn), (3) 
{2X* - l)iXn ~ X*) 



a;*(l — X*) 



To stabilize a fixed point by the ADP method, per- 
turbations Aa; are applied to the map's coordinate. The 
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value of the perturbation Ax is defined by the distance 
between the current system coordinate and the coordi- 
nate of the stabilized state: 



*^n+l ^*^n(l- '^n) ~t~ Ax,) 



(4) 



The ADP method is simple to use in practice. DiS'erent 
modifications of the adaptive method are therefore used 
in many papers devoted to experiments on the control of 
chaos. 

We consider two types of stabilization procedure: local 
stabilization and global stabilization. 

During local stabilization, the perturbations Ar and 
Aa; differ from zero only if the following condition is sat- 
isfied: 



< e. 



(5) 



Here e is a small value: we fixed e — 0.01. If the condition 
© is not satisfied then stabilization is absent, i.e. Ar — 
or Aa; = 0. 

During global stabilization perturbations are switched 
on when the condition (O is satisfied for the first time, 
and remain present for all future time. 

So, local or global stabilization involve modifications 
of the initial map (Q, and thus use another map in the 
form (PJ or Q). The fixed point a;* is an attractor of 
the new map. After the stabilization is switched on, a 
trajectory of the map tends to the fixed point a;*, and 
subsequently remains there. 

In the presence of noise the trajectory fluctuates in the 
vicinity of the stabilized state, i.e. noise-induced dynam- 
ics appears. In addition, noise can induce stabilization 
failures. For local stabilization they imply a breakdown 
in the condition (O , and for global stabilization they cor- 
respond to an escape of the trajectory from the basin of 
attraction of the fixed point x* . 

Our aim is to study these noise-induced stabilization 
failures and analyze the problem of how to suppress them. 
We therefore consider the maps Q and Q in the pres- 
ence of additive Gaussian fluctuations: 

Xn+i = (r + Ar„)x„(l - x„) + (6) 

. {2x* - l)ixn - X*) 

Ar„ = r r , 

a;*(l — X*) 

Xn+l = rXn{l - Xn) + AXn + DCn, (7) 
AXji — {x^ X ). 

Here D is the noise intensity; and ^„ is a Gaussian ran- 
dom process with zero-average (^) =0, delta-correlation 
function {^n^n+k) = S{k), and dispersion (^^) = 1. We 
use a high-speed noise generator 



II. CONTROL OF LARGE FLUCTUATIONS 

Large fluctuations manifest themselves as large devia- 
tions from the stable state of the system under the action 



of fluctuational forces. Large fluctuations play a key role 
in many phenomena, ranging from mutations in DNA to 
failures of electrical devices. In recent years significant 
progress has been achieved both in understanding the 
physical nature of large fluctuations and in developing 
approaches for describing them. The latter are based on 
the concept of optimal paths - the paths along which the 
system moves during large fluctuations. Large fluctua- 
tions are very rare events during which the system moves 
from the vicinity of a stable state to a state remote from 
it, at a distance significantly larger than the amplitude of 
the noise. Such deviations can correspond to a transition 
of the system to another state, or to an excursion along 
some trajectory away from the stable state and then back 
again. During such deviations the system is moved with 
overwhelming probability along the optimal path under 
the action of a specific (optimal) fluctuational force. The 
probability of motion along any other (non-optimal) path 
is exponentially smaller. In practice, therefore large fluc- 
tuations must of necessity occur along deterministic tra- 
jectories. The problem of controlling large fluctuations 
can thus be reduced to the task of controlling motion 
along a deterministic trajectory. Consequently, the con- 
trol problem can be solved through application of the 
control methods developed for deterministic systems 0| . 

Let us consider the control problem. Formally, the 
task that we face in controlling noise-induced large fluc- 
tuations consists of writing a functional R, the extrema 
of which correspond to optimal solutions of the con- 
trol problem, i.e. solutions with minimal required energy 
in 111 El- The form of the functional R depends on a 
number of different additional conditions related to e.g. 
the system dynamics, the energy of the control force, or 
the time during which it is applied 0, 0, ^3 ■ We will 
follow the work and consider the control of large fluc- 
tuations by a weak additive deterministic control force. 
Weakness means here that the energy of the control force 
is comparable with the energy (dispersion) of the fluctua- 
tions (see for details). In this case, the extremal value 
of the functional R for optimal control, which moves the 
system from an initial state a;' to a target state x-^ , takes 
the form Hal: 



i?opt(a;^F) = 5(°)(a;0±A5, 



(8) 



AS- = (2F) 



1/2 



N, 



k=Ni 



1/2 



where is the optimal fluctuational force that induces 
the transition from a;* to x^ in the absence of the control 
force; S*'"' is an energy of the transition, Ni and Nf are 
the tirnes at which the fluctuational force S,'^^^ starts and 
stops 24], and F is a parameter defining the energy of 
the control force. 

The optimal control force for the given functional 
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is defined by: 



.opt 



k=Nf 



1/2 



X d{Xn ^ ) , 



(9) 



where Xn"^"^^ is the optimal fluctuational path in the ab- 
sence of the control force. The minus sign in the ex- 
pression ^ decreases the probability of a transition to 
the state x-^ , and the plus sign increases the probabil- 
ity. It can be seen that the optimal control force 
is completely defined by the optimal fiuctuational force 
^^P', and the optimal fluctuational path Xn^°^^ , corre- 
sponds to the large fluctuation. Therefore to solve the 
control problem it is necessary, first, to determine the op- 



timal path X 



(O)opt 



leading from the state to the state 



c-^ under the action of the optimal fluctuational force 
^opt^ Thus, a solution of the control problem depends on 
the existence of an optimal path: it is obvious that the 
approach described should be straightforward to apply, 
provided that the optimal path exists and is unique. 

We consider below an application of the approach 
described to suppress large fluctuations in the one- 
dimensional map. The large fluctuations in question are 
considered here to correspond to failures in the stabiliza- 
tion of an unstable orbit. 

The control procedure depends on the determination 
of the optimal path and optimal fluctuational force and, 
to deflne them, we will use two different methods. The 
first is based on an analysis of the prehistory probability 
distribution (PPD) and the second one consists of solving 
a boundary problem for an extended map which defines 
fluctuational trajectories. 

The PPD was introduced in to analyze optimal 
paths experimentally in flow systems. We will use the dis- 
tribution to analyze fluctuational paths in maps. Note, 
that in 0, it was shown that analysis of the PPD al- 
lows one to determine both the optimal path and the op- 
timal fluctuational force. The essence of this first method 
consists of a determination of the fluctuational trajec- 
tories corresponding to large fluctuations for extremely 
small (but finite) noise intensity, followed by a statistical 
analysis of the trajectories. In this experimental method 
the behaviour of the dynamical variables Xn and of the 
random force ^„ are tracked continuously until the sys- 
tem makes its transition from an initial state to a 
small vicinity of the target state x-^ . Escape trajecto- 
ries x'^'^ reaching this state, and the corresponding noise 
realizations of the same duration, are then stored. 
The system is then reset to the initial state and the 
procedure is repeated. Thus, an ensemble of trajectories 
is collected and then the fluctuational PPD p!^ is con- 
structed for the time interval during which the system 
is monitored. This distribution contains all information 
about the temporal evolution of the system immediately 
before the trajectory arrives at the final state x^ . The 
existence of an optimal escape path is diagnosed by the 



form of the PPD p'^: if there is an optimal escape tra- 
jectory, then the distribution p'^ at a given time n has 
a sharp peak at optimal trajectory x°p*. Therefore, to 
find an optimal path it is necessary to build the PPD 
and, for each moment of time n, to check for the pres- 
ence of a distinct narrow peak in the PPD. The width 
of the peak defines the dispersion cr'^ of the distribution 
and it has to be of the order of the mean-square noise 
amplitude yHj The optimal fluctuational force that 
moves the system trajectory along the optimal path can 
be estimated by averaging the corresponding noise real- 
izations over the ensemble. Note, that investigations 
of the fluctuational prehistory also allows us to determine 
the range of system parameters for which optimal paths 
exist. 

To determine the optimal path and force by means of 
the second method we analyze extended maps 6, 7J using 
the principle of least action . Such extended maps are 
analogous to the Hamilton- Jacobi equation in the the- 
ory of large fluctuations for flow systems. For the one- 
dimensional map Xn+i = f{Xn) + D^n, the corresponding 
extended map in the zero-noise limit takes the form: 



Xn+l 
Vn+l 

9{Xn) 



fixn) + y„/g{xn), 

dXn 



(10) 



The map is area-preserving, and it deflnes the dynam- 
ics of the noise-free map Xn+i ~ f{xn), if Vn ~ 0. If 
yn 7^ then the coordinate Xn corresponds to a fluctu- 
ational path, and the coordinate y„ to a fluctuational 
force. Stable and unstable states of the initial map be- 
come saddle states of the extended map. So, the fixed 
point X* of the ADP O and OGY © maps becomes a 
saddle point of the corresponding extended map. Fluctu- 
ational trajectories (including the optimal one) starting 
from X* belong to unstable manifolds of the fixed point 
(x*, 0) of the extended map. 

The procedure for determination of the optimal paths 
consists of solving the boundary problem for the extended 
map (jTU)) : 



X — oci - - X 



= 
= 0; 



(11) 
(12) 



where x* is the initial state and x^ is a target state. 

To solve the boundary problem different methods can 
be used. For the one-dimensional maps under consid- 
eration, a simple shooting method is enough p^ . We 
choose an initial perturbation I along the linearized un- 
stable manifolds in a vicinity of the point (x* , 0) of the 
map (|10|) . The procedure to determine a solution can be 
as follows: looking over all possible values /, we deter- 
mine a trajectory which tends to the point (x^, 0). Note 
that, because these maps are irreversible there exist, in 
general, an infinite number of solutions of the boundary 
problem. The optimal trajectory (path) has minimal ac- 
tion (energy) S = Yl^=-oo Vn: here ?/„ is calculated along 
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the trajectory, corresponding a solution of the boundary 
task. 



III. NOISE-INDUCED FAILURES IN LOCAL 
STABILIZATION 

A breakdown of the condition jsj corresponds to a fail- 
ure of local stabilization, i.e. to the noise-induced escape 
of the trajectory from an e- vicinity of the fixed point x* . 
The target state x-^ corresponds to the boundaries of the 
stabilization region: x-^ — x* zt e. 

Instead of analyzing the maps ^ and {Tj) in the e- 
vicinity of the fixed point x* we can investigate linearized 
maps of the following form: 

Xn+i = ax„ + D^n, (13) 

here a is a value of derivative df{xn)/dxn in the fixed 
point X* . For the map © the derivative is equal to zero 
flOGY — 0, and for the map O cadp — —0.8. 

Let us investigate stabilization failure by considering 
the most probable (optimal) fluctuational paths, which 
lead from the point x* to boimdaries x*±e. For linearized 
maps the extended map can be reduced to the 
form: 

Xn+l = axn + —, (14) 
a 



with the initial condition (xq — x* , yo — 0) and the final 
condition x-^ = x*zte. It can be seen that a solution of the 
map H14|l increases proportionally to ?/„ = const/a" [l7l |. 
This means that, for the ADP map (O, the amplitude of 
the fluctuational force increases slowly but that, for the 
OGY map ©, the failure arises as the result of only one 
fluctuation (iteration). Because equation (fTHl is linear, 
the boundary problem will have a unique solution 'l6]. 
Thus, analysis of the linearized extended map H14I) shows 
that there is an optimal path, and it gives a qualitative 
picture of exit through the boundary x* ± e. 

Let us check the existence of the optimal paths through 
an analysis of the prehistory of fluctuations. To obtain 
exit trajectories and noise realizations we use the follow- 
ing procedure. At the initial moment of time, a trajectory 
of the map is located at point x*. The subsequent be- 
haviour of the trajectory is monitored until the moment 
at which it exits from the e-region of the point x* . The 
relevant part of the trajectory, just before and after its 
exit, are stored. The time at which the exit occurs is set 
to zero. Thus ensembles of exit trajectories and of the 
corresponding noise realizations are collected and PPDs 
are built. 

To start with, we will discuss these ideas in the context 
of the ADP map. Fig. 1(a) shows PPDs of the escape 
trajectories of the ADP map, and the corresponding noise 
realizations for the exit through the boundary (x* — e) 
are shown in Fig. 1(b). The picture of exit through the 



other boundary [x* + e) is symmetrical, so we present 
results for one boundary only. It is evident (Fig. 1) that 
there is the only one exit path. Note, that the path to the 
boundary (x* — e) is approximately 2.8x more probable 
than the path to the boundary (x* -I- e) . This difference 
arises from an asymmetry of the map in respect of the 
boundaries. 

Because for each boundary there is the only one exit 
path, the optimal path and the optimal fluctuational 
force can be determined by simple averaging of escape 
trajectories and noise realizations respectively. In Fig. 2 
the optimal exit paths and the optimal fluctuational 
forces are shown for the boundaries (x* — e) and (x* -|- e) . 
The paths and the forces coincide with a solution of the 
boundary problem (circles in the Fig. 2) of the extended 
linear map (|14|l . The time dependence of the dispersion 
all ol PPDs for the exit trajectories and noise realizations 
are shown in Fig. 3. As can be seen (Fig. 2) the optimal 
path is long, and the amplitude of the fluctuational force 
increases slowly, in agreement with analysis of the lin- 
earized map (|14|l . The dispersion cr^ of both trajectories 
and noise realizations decreases by construction as the 
boundary is achieved (Fig. 3). 

The optimal fluctuational force obtained (Fig. 2(b)) 
must correspond |l5j to the energy-optimal determinis- 
tic force that induced the stabilization failure. We have 
checked this prediction and found that the optimal force 
induces the exit from an e-region of the point x* : we se- 
lected an initial condition at the point x* and included 
the optimal fluctuational force additively; as a result we 
observed the stabilization failure. If we decrease the am- 
plitude of the force by 5-10%, then the failure does not 
occur. It appears, therefore, the deduced force allows us 
to induce the stabilization failure with minimal energy 
(see 10 for details). 

Using the o ptim al path and the force we can solve the 
opposite task |l3,[i3| — to decrease the probability of the 
stabilization failures. Indeed, if during the motion along 
the optimal path we will apply a control force with the 
same amplitude but with the opposite sign as the optimal 
fluctuational force has, then, obviously, the failure will 
not occur. Because we know the optimal force then, in 
accordance with the algorithm described above, it is 
necessary to determine the time moment when system is 
moving along the optimal path. For the ADP method the 
optimal path is long enough to identify that a trajectory 
is moving along the optimal path, and then to apply a 
control force. 

In the presence of control the map Q is modifled: 

x„+i = rx„(l - x„) + Ax„ + L»^„ + w„, (15) 
^Xji — (x^i X ) , 

here m„ is the deterministic control force. 

We use the following scheme to suppress the stabi- 
lization failures. Initially the control force is equal to 
zero {un = 0) and the map is located in the point x*; 
we continuously monitor a trajectory of the map H15|) 
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and define the time moment when the system starts mo- 
tion along the optimal path (a;„). We assume that the 
system moves along the optimal path (a;„) if it passes 
within a small vicinity of the coordinate {x-2) and then 
within a small vicinity of (x-a) (see arrows in Fig. 2(a)). 
Then on the following iteration we add the control force 
Un = -sign(^„)(^„), n = -1 (see Fig. 2(b)). 

In Fig. 4(a) dependences of the mean time (r) between 
the failures on the noise intensity D are plotted in the 
absence, and in the presence, of the control procedure. It 
is clear that the mean time (r) is substantially increased 
by the addition of the control, i.e. stability in the face of 
fluctuations is significantly improved by the addition of 
the control scheme. The efficiency of the control proce- 
dure depends exponentially jl2i] on the amplitude of the 
control force (Fig. 4(b)), and there is an optimal value of 
the control force, which is very close to the value (arrow 
in Fig. 4(b)) of the optimal fluctuational force. 

Now consider noise-induced stabilization failures for 
OGY map ©. An analysis of the linearized map has 
shown that the failure occurs as the result of a single fluc- 
tuation. We have checked the conclusion by an analysis 
of the fluctuational trajectories of the map much as 
we did for the ADP map. The optimal path and optimal 
force are shown in Fig. 5 for both boundaries, [x* + e) 
and (x* — e). An exit occurs during one iteration and 
there is no a prehistory before this iteration. It means 
that we cannot determine the moment at which the large 
fluctuation starts and, consequently, that we cannot con- 
trol the stabilization failures. The existence of a long 
prehistory is thus a key requirement in the control the 
large fluctuations. 

We can of course decrease the probability of a failure 
by increasing the e- region of stabilization. The maximum 
possible increase would correspond to inflnite boundaries 
- in which case we would be dealing with global stabi- 
lization. 



IV. NOISE-INDUCED FAILURES OF GLOBAL 
STABILIZATION 

To investigate fluctuational dynamics in the global sta- 
bilization regime, we consider the dynamics of the maps 
© and iQ with initial conditions at the fixed point 
xq = X* . We will first consider them in the absence 
of noise. The maps are shown on the plane a;„ — Xn+i in 
the Fig. 6. 

The map ((HJl (Fig. 6(a)) has three fixed points of pe- 
riod one: the point x* « 0.7368 is stable with the mul- 
tiplier /I = 0; the points X2 — Q and x\ w 0.5906 are 
unstable with multipliers fi « —3.04 and fi « 1.8016, re- 
spectively. The map has two attractors: the point x* and 
the attractor at infinity [2^. Their basins of attraction 
(Fig. 6 (a)) are self-similar (fractal) The point 

x\ and its pre-images by backward iteration lie on the 
basin boundaries of the attractors |23 . In the intervals 
X e (-0.183,0.5906) and x e (0.862, 1.027) the basins of 



the attractors alternate and are of different length. The 
interval x G (0.5906, 0.862) corresponds to the widest 
basin of the fixed point x* . The boundaries of this basin 
are defined by the unstable point x^ and its pre-image 
x{* . The semi-infinite intervals x G (—00,— 0.183) and 
X e (1.027,00) correspond to basins of the attractor at 
infinity. The boundaries of the semi-infinite intervals are 
defined by the points x~°° = -0.183 and a;°° = 1.027, 
which correspond to the cycle of period 2. 

The map Q (Fig. 6(b)) has two fixed points: the point 
X* « 0.7368 is stable with multiplier fi w —0.8; and the 
point w 0.2632 is unstable with multiplier /i w 2.8. 
The map has two attractors: the fixed point x* and the 
attractor at infinity. The basins of attraction are smooth 
(Fig. 6 (b)). The first boundary of the basins is the point 
x* and the second boundary is a pre-image x{* of the 
point xl- 

So, each of the maps has two attractors, but the struc- 
ture of their basins of attraction are qualitatively differ- 
ent. 

We now consider these maps © and ^ in the presence 
of noise. Noise can induce escape from the basin of the 
fixed point x* , corresponding to failure of the stabiliza- 
tion. As before we examine the dynamics of the escape 
trajectories obtained for extremely small noise intensity 
in order to determine the optimal path and the optimal 
force. Fluctuational escape trajectories of the map {Tj) 
are shown by dots on the plane [xn — Xn+i) in Fig. 6(b). 
As can be seen, there is one escape path, and the es- 
cape trajectories pass through the unstable point x^. In 
Fig. 7 the optimal path and the optimal force obtained 
by averaging the escape trajectories and noise realiza- 
tions respectively are shown by crosses. The stabiliza- 
tion failure clearly possesses a long prehistory. From the 
point of view of the control procedure, the presence of a 
large deviation of the system coordinate (xn) at the time 
moment n = — 1, and the smaller deviation of the fluctu- 
ational force (^„) at the next time moment (n = 0), are 
important. This is because the first fluctuation of coordi- 
nate Xn can easily be identifled and distinguished from 
non-optimal fluctuations in the vicinity of the stable state 

X*. 

Next, we examine the process of escape for the map 
©• Fig. 8(a) shows escape trajectories superposed at 
the time moment when the trajectory crosses the basin 
boundary at the point x~°°. It is evident that there 
is no selected escape path. The escape trajectories can 
be divided into several groups with different probabili- 
ties. With maximum probability (almost 50%) the escape 
trajectories follow the arrowed path in Fig. 8(a)) corre- 
sponding to motion in the direction of the point x~°° 
without any jumps in the opposite direction. The other 
paths include jumps in the opposite direction. The width 
of the distribution of fluctuational paths is comparable 
with the noise amplitude and there is no a specific fluc- 
tuational force. In Fig. 6(a) the escape trajectories are 
shown on the plane {Xn - Xn+i)- Is can be seen that, af- 
ter the point xl, the escape trajectories are located close 
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to the trajectories of deterministic map, so we can sup- 
pose that after the point the motion has the character 
of directed diffusion. The interval between the points xl 
and Hes within the fractal basin, and this fact im- 
plies a variety of escape paths. Indeed, within a small 
vicinity of the point xl there is a piece of basin of the 
attractor at infinity. For escape, therefore, it is enough 
to bring the trajectory only to this basin. However, the 
size of this basin is small and a weak fluctuation can of 
course move the trajectory back to the basin of point x* 
and vice versa. As a result, the trajectory can spend a 
long time in the vicinity of the point xl: it can return 
to the point x* , as well as escape from the basin of the 
point X*. 

Thus, the fractal structure in the basin of attraction 
leads to complex behaviour of the escape trajectories; 
they can spend a long time in the fractal basin; motion 
in the direction of the attractor at infinity has the largest 
probability. 

Investigations of escape from the point x* to the vicin- 
ity of the point xl have shown (Fig. 8(b)), that there is 
no specific path within this interval, so that we cannot 
determine the optimal path or the optimal force using 
an analysis of the escape trajectories. It is possible to 
select several different favoured paths (thick lines in the 
Fig. 8(b)), but dispersion of the trajectories for each of 
them is much larger than the noise intensity used. 

We now determine the escape optimal paths and the 
forces by solving the boundary problems Hll|l and 112|l 
for the extended maps: 



and 



Xn+l 
yn+1 
f{Xn) 



Xn+l 
Vn+1 
f{Xn) 

9{Xn) 



fiXn) + Vn/giXn), 
Vn/giXn), 

{r + Ar„)a;„(l - Xn), 

dXn 



f{Xn) +yn/g{Xn), 

yn/g{xn), 

rXnil - Xn) + AXn, 
dfjXn) 
dXn 



(16) 



(17) 



which correspond to the maps lO and iQ. In such a 
way we have used the extended map H14() to analyze the 
linearized map l(T^ . 

First, we consider the results of solving the boundary 
problem for the extended map H16|) . To do so, we use 
a shooting method, with boundary conditions Hll|l and 
(fT^ . where x^ = x*. Since the derivative g{xn) = ^'^g^"^ 
of the map © at the point x* is equal to zero, we can- 
not calculate eigenvectors of the point (a;*, 0) of the map 
(|16|l . Therefore, as a parameter of the boundary prob- 
lem we choose an initial perturbation yg, since it defines 
all the trajectories going away from the point (a;*,0). 



Four solutions of the boundary problem, obtained nu- 
merically, are found to have practically the same action 
S. Four escape paths and noise realizations {ti — t^) of 
the map l|16|) corresponding to these solutions are shown 
in Fig. 9. The trajectory has the minimum activation 
energy S « 0.0115 and the energies of other trajectories 
are practically the same: S « 0.0123. All the optimal 
trajectories lie on a stable manifold of the point {xl,0), 
and the stable manifold goes to the point {x* , 0) (Fig. 10). 
If we take into account the fact that the noise intensity is 
finite during the experimentally analysis of escape trajec- 
tories (Fig. 8), then the fluctuational trajectories of the 
map (jSJ form a wide bunch around the optimal paths 
and trajectories can go along different the optimal paths 
at different time intervals. Thus, for the OGY map jnj, 
the only way to determine optimal paths and forces is by 
solution of the boundary problem for the extended map, 
whereas an analysis of the PPD is not successful. 

Now, let us consider the solution of the boundary prob- 
lem for the map H17(l . We have defined an unstable direc- 
tion of the point (a;*,0) and used the length of a vector 
I along this direction as a parameter of the boundary 
problem. There is just one solution for which the value 
of action S = 0.0449, which is slightly smaller than the 
value S — 0.493 calculated by using the PPD. The cor- 
responding optimal path and optimal force are shown in 
Fig. 7 together with the path and the force found by us- 
ing PPD. It can be seen (Fig. 7), that the optimal paths 
and the forces obtained by calculated PPD and by using 
the extended map are practically the same. 

Thus, we have defined the optimal path and the op- 
timal force corresponding to global stabilization failures, 
and we have compared two methods for determination of 
the optimal path and force: the first method being based 
on an experimental analysis of the prehistory probability 
distribution, and the second one being based on solving 
the boundary problem for an extended area-preserving 
map. The latter method allows us to determine the opti- 
mal path and force for both the maps © and {Tj) whereas 
experimental analysis of prehistory probability is only 
successful for ADP map (|7|). 

Because there is no an unique escape path for the OGY 
map, it is impossible to apply the algorithm described 
above for controlling stabilization failures. We note how- 
ever that, since we know the dynamics of the fiuctua- 
tional trajectories, it is still possible to realize control of 
the fiuctuations by using another approach. For example, 
a control force can be added whenever the system comes 
to the vicinity of the point xl. In this case, however, the 
size of the vicinity and the magnitude and form of the 
control force are ill defined. 

For stabilization of the ADP map, the opposite situ- 
ation applies: there exist an unique optimal path and a 
corresponding optimal force. Consequently, we can real- 
ize a procedure for the control of large fluctuations. It 
is similar to that described above for local control. We 
monitor trajectories of the map (0 to identify the large 
deviation ((a;„), n = —1 in Fig. 7) and in the next it- 



7 



eration we add the control force u„ = — (^n), n = 0. 
The dependences on noise intensity D of the mean time 
(r) between stabiHzation failures in the absence and in 
the presence of control are shown in Fig. 12(a). The de- 
pendence of (r) on the amplitude of the control force is 
shown in Fig. 12(b). The suggested control procedure is 
evidently effective. 

Conclusion 

We have considered noise-induced failures in the sta- 
bilization of an unstable orbit, and the problem of how 
to control such failures. In our investigations, they cor- 
respond to large deviations from stable points. We have 
examined two types of stabilization, local and global, and 
therefore analyzed fluctuational deviations of different 
size. We have shown that, for local stabilization, noise- 
induced failures can be analyzed effectively in terms of 
linearized noisy maps. 

Large noise-induced deviations from the fixed point 
in one-dimensional maps have been analyzed within the 
framework of the theory of large fluctuations. The key 
point of our consideration is that the dynamics of the 
optimal path, and the optimal fluctuational force, corre- 
spond directly to stabilization failures. We have applied 
two approaches - experimental analysis of the prehistory 
probability distribution and the solution of the bound- 
ary problem for extended maps - to determine the op- 
timal path and the optimal fluctuational force, and we 
have compared their results. For local stabilization, the 
two approaches give the same results. For global stabi- 
lization, however, the solution of the boundary problem 
enabled the optimal path and optimal fluctuational force 
to be determined for both the OGY and ADP maps, 
whereas investigation of fluctuations' prehistory gave the 
optimal path and force for the ADP map only. 

A procedure for the control of large fluctuations in one- 
dimensional maps has been demonstrated. It is based on 
the control concept developed in 12] for continuous sys- 
tems. We have introduced an additional control scheme 
which significantly improves the stabilization of an un- 
stable orbit in the presence of noise. It was successful 
for the ADP method of stabilization, and problematic 
for the OGY method. We have shown that the control 
procedure has limitations connected with the existence of 
unique optimal path and the presence of long time prehis- 
tory of large fluctuation. The relationships between the 
large fluctuation dynamics and the control procedures are 
summarized in Table I. 

The considering of the control problem is relevant to 
a continuous system which has one-dimensional curve in 
the Poincare section, for example, to the Rossler system. 



For such a continuous system we can formulate the task of 
control as a control in discrete time moments (moments 
of intersection of the Poincare section) by using impulse 
actions. Intervals between the moments were used to 
calculate and to form a control force. Note, that the 
similar approach is wide used in control technique. 

The present control approach has limitations which 
consist of necessity of studying the fluctuational dynam- 
ics of systems before considering of the control problem. 
Such studying can be done with use of extended maps of 
system, if it is known a system model, and/or experimen- 
tally by fluctuation prehistory analysing. For the local 
stabilization a system model can be easily written down 
by determination of eigen-value of a stabilized unstable 
point; there are a lot of effective methods to do it [2ll |. 
For the global stabilization there is no such a way and 
we need to investigate the fluctuation prehistory. Our 
investigations have shown that in this case we can met 
problems of determination of the control force. Indeed, 
we have shown that for the global stabilization of OGY 
map there are several most-probable escape paths with, 
practically, the same energy. As the result, a real escape 
path can be a combination of the different most-probable 
paths, so an escape trajectory does not go along a de- 
fined path as for the ADP map. Moreover, we can not 
determine the fluctuational force and, consequently, the 
control force, since we use for it averaging fluctuational 
trajectories which follow along an unique path. So, the 
control procedure is inapplicable. It is obviously, that for 
succesful control we have to change the control strategy. 
For example, we can try to predict a fluctuational action 
locally whereas now we try to know the full fluctuational 
dynamics. The local prediction can be based on a combi- 
nation of real time prehistory analysis and reconstruction 
of the extended system 22] . 

Additionally, noise-induced escape through fractal 
boundaries has been studied in a one-dimensional map. 
It was found that fluctuational motion across fractal 
basins has a non-activation character. It was also estab- 
lished that there are several optimal escape paths from 
the fixed point of the OGY map © whereas, for the ADP 
map ([7|l, the escape path is unique. We infer that the ex- 
istence of several paths in the OGY map © is connected 
with the fact that the stable manifold of the boundary 
point (a;i,0) goes to the fixed point (a;*,0). 
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FIG. 1: PPDs of the exit trajectories (a) and noise re- 
alizations (b) of ADP map for the boundary [x* — e). The 
thick dashed hnes indicate e-region of stabilization. The thin 
dashed lines connect maxima of PPDs. The noise intensity is 
D = 0.0011. 




FIG. 2: The optimal paths (a) and the optimal forces (b) for 
exit through the boundary {x* —e) (solid line) and the bound- 
ary [x* -\- e) (dashed line) for ADP map. Circles indicate the 
optimal paths and forces obtained by solving the boundary 
problem for the linearized extended map dJ. The optimal 
paths and forces used in the control procedure are marked by 



arrows. 
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FIG. 4: (a) The dependences of mean time (r) between sta- 
bilization failures on noise intensity D in the absence (circles) 
and in the presence (crosses) of the control. The size of the 
stabilization region is e = 0.01. (b) The dependence of the 
mean time (r) on the amplitude of the control force u„ is pre- 
sented for the ADP method. The value of (r) corresponding 
to the optimal fiuctuational force is marked by the arrow. 




FIG. 5: For the OGY map, the optimal path (a) and the 
optimal force (b) arc shown for exit through the boundary 
(x* — e) (crosses) and the boundary {x* + e) (circles). 
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FIG. 6: The OGY map (a) and the ADP map (b) on the 

plane (.Xn — .x'n+i) arc shown by the thick solid line. Basins 
of attraction of the fixed point x* (white regions) and the at- 
tractor at infinity (black regions) arc shown at the bottoms 
of the figures. The dashed lines indicate locations of the fixed 
points of the maps and the points defining the basin bound- 
aries. Escape trajectories are shown by dots. The thin line in 
figure (b) corresponds to the optimal path. 




FIG. 7: The optimal path (a) and the optimal force (b) 
obtained by experimental analysis of the PPD for the ADP 
map Q (crosses) and by solving the boundary problem for 
the map (circles). 





FIG. 8: 100 escape trajectories of the OGY map collected in a 
vicinity of (a) the point x~°° and (b) the point x^. The size of 
the vicinity is defined as the mean square of noise intensity D. 
The dash-dot lines indicate the location of the boundary point 
xl and its pre-image x{*; the dashed lines in (a) represent 
boundaries of the basins with fractal structure, i.e. of the 
points x~°° and x°° . The thick lines in (b) correspond to 
different escape paths; the grey dots are coordinates of escape 
trajectories. The noise intensity is D — 0.018. 




FIG. 10: The stable manifold of the point Xj of the OGY 
map. The symbols indicate the different optimal paths, using 
the same coding as in Fig. 9. 
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FIG. 11: (a) Dependences of the mean time (r) between sta- 
bilization failures on noise intensity D in the absence (circles) 
and in the presence (crosses) of the control. The stabiliza- 
tion is global, using the ADP method, (b) Dependence of the 
mean time (r) on the amplitude of the control force u„. 



